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Abstract 

Combining three mechanisms, we reanalysis processes of B ^ rj' K{K*),'qK{K*) and 
calculate their branching ratios. The results are compared with other mechanisms in 
the literature. The striking feature of the gluon fusion mechanism is emphasized and 
its experimental test is discussed. 
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1 Introduction 



Since the CLEO collaboration reported large branching fractions |I|]: 

Bn{B^^ri'Xs) = (6.2 ± 1.5 ±1.3) X 10-^(2.0 < P^/ < 2.7GeV"), (1) 
Bn{B^ r]'K^) = {Q.5tli ± 0.9) x 10"^ . (2) 

It has received much attention because it is a great experimental achievement in rare B decays 
which are dominated by penguin contributions and the branching fractions are surprisingly 
large compared with earlier theoretical estimation. 

Now it seems that we have a qualitatively reasonable interpretation of the semi-inclusive 
decay B — > rj'Xs: It is due to the special property of t]' which has anomalously large coupling 
to gluons. But as to the exclusive decay B —>■ rj'K there are still several possible mechanisms. 
This is partly because it is less reliable to estimate exclusive decay than to semi-inclusive 
decay due to our complete ignorance of the hadronization process. 

In the following we first carefully estimate the contribution from the conventional mech- 
anism. We find that the terms, such as ( r^' | S75S | 0) ( K~ \ sh \ B~) , are dominant terms 
if using Dirac equation to estimate them, but the large factors such as have large uncer- 
tainties because the quark mass is uncertain, so we also use PQCD method to estimate them 
and get comparatively small result which indicates that we need new mechanisms to account 
for the experiment result. Secondly we discuss some possible new mechanisms. We find that 
if 6 ^ CCS —>■ ST]' is the dominant mechanism, the branching fraction of B ^ rj'K* would 
be too large to fit the CLEO data. Thirdly we try to combine the mechanisms proposed in 
ref [6] (two gluons fusion) , ref [7] (non-zero f^, contribution) and the conventional mechanism 
to compute the branching fractions of B ^ ri'K{K*),riK{K*). The results are compared 
with other mechanisms in the literature. 

The article is organized as follows: Section 2 gives the theoretical formula for the decay 
amplitude of 5 — > r]'K using factorization assumption and Dirac equation, and some uncer- 
tainties in the formula are briefly discussed. Section 3 is devoted to the PQCD method to 
estimate the hadronic matrix elements such as ( r^' | S75S | 0) ( \ sb \ B^) . In section 
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4, some possible mechanisms are discussed and used to fit the data. Section 5 is for the 
concluding remarks. 



2 Preview 

The standard theoretical frame to estimate the non-leptonic B decays is based on the effective 
Hamiltonian 



H, 



eff 



V2 



10 



1=3 



(3) 



( \i = VibV*.) and BSW model 0. But there exists large uncertainties in estimating the 
hadronic matrix elements like 



( 7]' I S75S I 0) , { K- \ sb \ B-) . 



(4) 



In the literature Dirac equation is used to estimate such matrix elements: 



( v' I s'j^s I 0) 
K- \sb\ B-) 



K- I 57^6 I B-) 



nib - ms 

With these relations one can write down the amplitude of rj'K^ |^: 



(04 + aio + 2(a6 + 



m 



K 



{rus + mu){mb - rriu) 



+ (2a3 - 2a5 - ^07 + ^ag)X2u + (03 - ag - ay + ag)X2c 

1 1 1 / In ^l' 

+ (03 + 04 - 05 + o«7 - o«9 - o«10 + (Oe - 7^«8) 7 ^ ^>X2s 

2 2 2 2 nislnib — iris) 



m 



+ (04 + oio + 2(a6 + 03)7 ■ r)X3 

{iris - mu){mb + mu) 



}• 



where X, are factorizable terms 



(5) 
(6) 



(7) 



Xi ={K-\ {su)v-A I 0) ( V I {ub)v-A I , 
^2, ={v'\ iqq)v~A \0) {K-\ {sb)v^A I B-) , 
X3 = ( r]'K' I (sM)y_A I 0) ( I I B-) 



and 



{ K- I s^,u \0){v'\ub\ B-) = — ^ -X,, 
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iK- I I 0) ( I I B-) = — -X3. 

[rris - mu){mb + m^) 



But the large factors such as in the amphtude have large uncertainties because the quark 
mass is uncertain. 

Through the numerical calculation, we find Equation(7) gives a branching ratio of 
3.5 X 10~ , and if we neglect the term ^ (m^-m in the equation(7), the result is 
6.63 X 10^^, which is too small. So we find that such a term is very important, and we 
try to find an alternative way to estimate such matrix elements. 



3 Perturbative QCD calculation 

Brodsky et al. [^] has pointed out that because of the large momentum transfers involved in 
the decays of B to light mesons, the factorization formula of PQCD for exclusive reactions 
becomes applicable: the amplitude can be written as a convolution of a hard-scattering 
quark-gluon amplitude T/j, and meson distribution amplitudes Q"^) which describe the 
fractional longitudinal momentum distribution of the quark and antiquark in each meson. 
An important feature of this formalism is that, at high momentum tranfer, long-range final 
state interactons between the outgoing hadrons can be neglected. In the case of nonleptonic 
weak decays the mass of the heavy meson Mfj establishes the relevant momentum scale 
q2 ^ M^, so that for a sufficiently massive initial state the decay amplitude is of order 
even without including loop corrections to the weak hamiltonian. The dominant 
contribution is controlled by single gluon exchange. 

We try to use PQCD method to estimate those hadronic matrix elements such as ( 77' | 
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575* I 0) and { K \ sb \ B ) . The wave functions of B and K are: 

- j=,^,MyWK + , (10) 

where Ic is an identity in color space. In QCD, the integration of the distribution amphtude 
is related to the meson decay constant 

J (f)K{y)dy = ^^fK , J (f)B{x)dx = ■ (11) 

Then we can write down the amphtude of Fig.l as 

( K~ I S7^75W I 0) PQCD = 3 X J dy(j)K{y)Tr [75 (^k + Mxhi^-f^] = /kQh , 
( K- I S75« I 0) PQCD = 3 X / dycPK{y)TT [75 ((/^ + Mk^^,] = J^Mk • (12) 

Compared with the widely used results 

( K~ I S7^75?i I 0) = -i /fc^f^, 

(ms+mu) 



{ K I S75M I 0) Dirac = iu;S;r-)fK, (13) 



so, except for a — i factor, ( K | sjfj^j^u \ 0) pqcl* = { K \ sj^j^u \ 0) , but ( K \ sj^u \ 
0) PQCD { K~ \ s^5U I 0) Dirac ■ Similarly, we can get 

( V' I Q7^l5q I 0) PQCD = ( I 575? I 0) PQCD = f^'Mrj>. (14) 

In a consistent way, we can use perturbative QCD to estimate the matrix elements like 
( K~ I sjfj,b I B~) and ( K~ \ sb \ B~) (Fig.2, Fig. 3), we have neglected the fermi motion 
of quarks, while the gluons in the Fig. 2, 3 are hard because 

= {xPb - (1 - y)PKf ^ - y)Ml ~ IGeV (15) 

so, we can use perturbative QCD method to calculate the amplitude and it turns out to be 

2 f 

{ K~ I s-f^b I B~) PQCD = -gfl' y dxdy(f)B{x)(f)K{y) (16) 
^ k'^P^ k%P^ -ml) ^ ' 



In order to get the quantitative estimation, we take the wave functions as 

^b{x) = 0^5{x - ee) , ^K{y) = \j^fKy{l - y) ■ (17) 

we get 

( K I 57^6 I B ) PQCD = Y^^JBJK I ^2^4 + 

[ , -4|/M| + SyMeMK + Sm^MB - AMbttik \ , , , 

/ ^yy 2 /i^2/ — 2^ ip + 

7 -t%M'B{yMf^-nii) J 

= nP + Q)f^, (18) 

where we have neglected another term F' {p — q)^ which has neghgible contribution to the 
amphtude ( rj' \ 57^^755 | 0)( \ s'-ff^b \ B^). For simphcity, we will drop such terms like 
F' [p — q)^ in the following, since their contributions always proportion to m^,, m\ etc. The 
numerical results for F are presented in Table 1, where we have taken 

a,((A;2)) - 0.38, = 200Me\/, fx = IGOMeV. 

From Table 1, we find that our PQCD results are sensitive to the values of parameter 
and rrih, and seem small compared with the BSW result: 

( K- I s^^b \B-) = Fi^{ml,){p + q), = 0.35(p + q)^. (19) 

In fact the PQCD results are comparatively small in many cases, however, we believe that 
the ratio 





S75S 


1 0) pqcd{ K- 


1 sb 


B ) PQCD 






0) pqcd{ K~ I 


s'~)^^b 


B ) PQCD 



is more reliable because main uncertainties can be canceled. We can write down the matrix 
element ( K'^ \ sb \ B^) as 

{ K \ sb \ B ) PQCD = --9 J dxdy(j)B{x)(j)K{y){ 

Tr [75(9/+ MK)(ff'b + mh)r(ff' + Mb^^Iu] , 

Stt ^ ^ ,MkMUI - 2eB) - Mb{MI, - fM|) 



^--sj^jnv ^2^^4 

f ArubMBMK - m^Ml - (1 + y)MKMl + 8M^ ^ 
^1 -eBMl{yMl-ml) ^'^^ 



With eq.(14), (18), (21), we find that the ratio is really not sensitive to the values of param- 
eter eb and rri},. 

i?m = (0.19 ~ 0.20)(0.05 <eB< 0.08, A.SGeV <mi, <^.QGeV) (22) 
If our results are right, the ratio where Dirac equation is used, 

0.9 (23) 



{i\ 




0) Dirac{ K 


1 1 


^ ) Dirac 


{ri'\ 




0) ( i^- 1 




\B-) 



2ms(mb - rus) 
{ms{mb) ~ lObmev, 4:.8GeV < rub < b.OGeV) 

may be misleadingly large. 

It is difficult to quantitatively estimate the hadronic matrix element { rj' \ ub \ B~) 
because we don't know the wave function of r]' . Fortunately in the calculation of ( K~ \ sh \ 
B~) we find that ( K~ \ sb \ B~) pqcd ~ ( K~ \ sb \ B~) Dirac , so we can reasonably 
assume that { rj' \ ub \ B~) pqcd ^ { rj' \ ub \ B~) Dirac and get the estimation: 







S75U 


0) pqcd{ v' I 


ub 


B ) pqcd 


{K-\ 


s^i^^^u 


0) pqcd{ v' 


1 u-ff,b 1 


B ) PQCD 



0.1. (24) 



Because the term 



( K- I -s^,u I 0) ( r?' I ub I B-) = — — ^ -Xr (25) 

{ms + mu)[mb - mu) 

in eq. (7) (using Dirac equation) only contributes about 10 % to the total amplitude, our 
PQCD estimation (eq.24) will not substantially change the total amplitude though it is not 
very precise. 

For the hadronic matrix element 

2 

( rj'K~ \ su \ 0) (01 ■U756 I B~) — —- -X3 (using Dirac equation) (26) 

{ms - mu){mb + mu) 

It contribute about 10 % to the total amplitude in eq.(7) although it has an uncomfortable 
large factor . "^^ s , so we can still neglect it safely. 

From the above analysis, we obtain the revised conventional amplitude of B~ — > r]'K~ 

(assume f^, = 0) 

Mr.ciB- ^V'K-) = ^{VubV:^{affX, + al^^X,u) 



-VtbVZ [(04 + aio + 0.2 (ae + a^))X^ 

1 1 
+ (2a3 - 2a5 - -a^ + -agjAa^ 

+ (03 + 04 - as + - ^ag - ^aio + 0.20(2a6 - a^))X2s 



}(27) 



Taking VthV^ = -0.039, we get 

Bn{B' r]'K-) ~ 0.9 x 10"^ (28) 

(which is much smaller than the result of 3.5 x 10^^ where Dirac equation is employed.) So 
we find that the conventional decay mechanism can not account for the experimental data. 
We need some new mechanisms to interpret the experiment. 

4 Some new mechanisms 

There are several mechanisms which may enhance the decay rate of B"^ v'^^^ they all 
use the unique properties of 17' : 



(i) b — s> sg* and g* gt]' ^ via QCD anomaly , it is important to the inclusive decay 

B i]'-^s, but it seems difficult to realize its contribution to two— body exclusive 
decay B^ rj'K^. 

(ii) b sgg srj' [Q. Because the effective vertex oib sgg is very small, it is impossible 

to account for the large B^ t]'K^ branching ratio. 

(iii) Zhitnitsky have proposed a hopeful mechanism that rj' can be directly produced 
through b ccs — > rj's if ( r^' | cq^'~^^c \ 0) = —if^,P,l^, 7^ . Unfortunately we can't 
determine f^, precisely from QCD. If this mechanism dominates over other mechanisms 
of the decay B^ rj'K^ , Zhitnitsky gets f^, = lAOMeV. However one would find 

(29) 



Bn{B- 7]'K-) 



X' 



2c 



2c 



0.83. 
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and the branching ratio for B"^ rj'K*^ would exceed the upper hmit of the experi- 
mental data. Furthermore, if we take the contribution from the conventional mecha- 
nism into account we find that the branching ratio BTZ{B^ —>■ rj'K*^) is larger than 
BTZ{B^ rj'K^), this strongly disagrees with the experiment [|l|. So it is unlikely 
to be the dominant decay mechanism. In the following we tend to treat f^, as a free 
parameter in the range of —40 MeV to 40 MeV to take into account the contributions 
of this mechanism. 

(iv) Another possible mechanism is proposed in 0. This mechanism is motivated by the 
fact that both the recoil between rj' and K" and the energy released in the process 
are large. The gluon from b sg vertex would carry energy about Mb/2 and then 
materializes to t]' and emits another hard gluon to balance color and momentum. This 
is a hard scattering process, and the perturbative QCD method is applicable. This 
mechanism is depicted in Fig. 4 . We re-examine its contributions and the numerical 
results are presented in Table 2. And we find that this mechanism is very important. 

We do not discuss some other new mechanisms which involve new physics beyong the 
SM. We think that the contributions of the SM should be carefully examined first. So, in 
what follows, we try to combine the mechanisms (iii) and (iv) with the revised conventional 
amplitude eq.(28). The amplitude of mechanism (iv) is P 

M^v = ^alC,ffCp32{V,,V,l) 



, F^kflp ■ kiq ■ k2 - p ■ k2q ■ h] + F:^MBmb[q ■ k2kj - q ■ kiki ■ 
X / dxdy(f)B{x)(j)K{y) _ fcafcifcl ' ^ ^ 



KA — 



and the amplitude of mechanism (iii) is 

V^V:,al^^X2c - VthVL («3 - as - ar + a^)X2} (31) 
The total amplitude is 

Mtotal = Mr.c + Mm + Miy (32) 

The numerical results are listed in Table 2, where we have taken /§ = 1.38/,r, /o = 1-06/7^. It 
is shown that a nonzero may be important to explain the experimental data. From Table 



2, one may note that f^, = —30MeV seems fit the experiment better than = —15MeV. 
However, considering the uncertainties in estimating the branching ratios, for example the 
PQCD method often gives comparatively small result, we choose = —15MeV as an input 
parameter to estimate some other related channels, such as 

B-^ri'K*-, B-^7]K-, B-^T]K*-, 



B" r]'K*\ B" r/ir^ B" r]K 



*o 



We present a complete expression of the amplitude of -B — > r]'K*^ as an example: 



M 



total 



(33) 



where 



M 



with 



IV V* (n'^ff Y' -U n'^ff Y' 



VtbV* 



[a4 + aio + 0.2(a6 + a8))X'i + (2a3 - 2a5 - + ^ag)X2^ 



+ (as + 04 - 0^5 + ^a? - ^ctg - ^ctio + 0.2(2a6 - a8))^2s 



M 



V2 



ycbV*^aJ^ X'^c - VtbV^Xo.^ - 05 - 07 + aQ)X. 



2c 



(34) 
(35) 



X[ = {K*-\ {su)v-A I 0) ( r/' I {uh)v^A I B') 



X. 



2q 



{ rl I mv-A I 0) ( K*- I {sh)y_A I B-) . 



In the spirit of |^, we can calculate We take the wave function of K* as 

1 l(j 

where e is polarization vector of K* and after a direct calculation, we get 
M',, = ^a2^e//C^32(l/,,l^,*) X j dxdyM^)<pK4y) 



. ^ M bKc ■ fci)(g ■ fca) - (e ■ A:2)(g • h)] + M^4(e ■ k2){p ■ h) - (e ■ h)ip ■ ^2)] 

fci ■ k2klk2 



} 



(36) 



(37) 



, ^ [(e • ki)ip ■ ki){q ■ k2) - (e ■ A;2)(p ■ ki){q ■ ki)] + MbMk* [-(e ■ ki){ki ■ ^2) + (e ■ ^2)^;?] , 
+^2^bi r — 7-7:27:2 )• 
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The numerical results are shown in Table 3, where we have taken fx* = ISOMeV and 
eB = 0.06(0.07). It is shown that all predictions on the decays are under the upper limit 
of experimental data. Compared with the results of and [|I^ ,we have very different 
predictions on B~ r]' K*~ and ■q'K*^: 

Bn{B- T]'K*~) = 3.18(2.61) x 10~^ Bn{B^ r]'K*^) = 3.33(2.74) x 10"^ (38) 



which is two orders larger than the predictions of and [|I0|- If we use conventional 
amplitude used in the literature and take into account the contribution of Mechanism 
(iv) 1^, the results are listed in the sixth column of Table 3. We find that the branching 
ratios BTZ{B^ — > rj'K^) and B1Z{B^ — > rj'K^) all agree with the CLEO results. However, 
the predictions for B^ rj'K*" and B^ —>■ rj'K*^ are 

BTZ{B- r]'K*-) = 7.41 x 10-^ Bn{B^ t]'K*^) = 7.29 x 10-^ (39) 

where = 0.06, which are still one order larger than other predictions |§|,|13]- This character 
due to the mechanism (iv) can be tested in the near future. 



5 Concluding Remarks 

In this paper, we first carefully examine the conventional decay mechanism (using standard 
effective weak Hamiltonian and BSW model) which contributes to the exclusive decay B^ 
ri'K~ . Instead of using Dirac equation, we use an alternative PQCD method to estimate the 
hadronic matrix elements like ( r/' | S75S | 0) ( K~ \ sh \ B~) and get comparatively small 
results. We find that the standard theorectical frame can only give Br{B^ —>■ r]' K~) ~ 
1 X 10^^ , which is too small to account for the experimental data. 

We also discussed some new mechanisms within the SM, we find that f^, = lAOMeV 
will predict a large branching ratio of 5° rj'K*^ , which is in disagreement with the 
experimental result, so the mechanism through h ccs srj' is unlikely to be dominant 
although its contributions to the decay B —>■ r]'K{K*) are important, we combine the 
mechanisms of ref [0], ref 0] and the conventional mechanism to interpret the experimental 
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data. We find that the numerical results of BTZ{B^ rj'K^) are in good agreement with 
the experimental results and the predictions on the other decay channels are all under the 
experimental upper limit. It is interesting to note that we predict large branching fractions: 

Bn{B- riK*-) = 3.18(2.61) x 10-^ 
i37^(5° ^ T]'K*'^) = 3.33(2.74) x 10"^ 

which are still under the experimental upper limit but much larger than the predictions 
of and |10[. This can be tested in the near future. From Table 3, we find that the large 
results mainly come from the two gluons fusion mechanism which predict the same order 
branching fraction for B —>■ rj'K and B — > rj'K* 
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rrib = b.OGeV 


ee = 0.05 


SB = 0.06 


es = 0.07 


SB = 0.08 


F 


0.31 


0.24 


0.20 


0.17 


Rm 


0.190 


0.192 


0.193 


0.195 




rub = A.9GeV 


es = 0.05 


SB = 0.06 


SB = 0.07 


es = 0.08 


F 


0.28 


0.22 


0.18 


0.16 


Rm 


0.192 


0.194 


0.196 


0.198 




rrib = A.SGeV 


es = 0.05 


SB = 0.06 


eB = 0.07 


SB = 0.08 


F 


0.26 


0.20 


0.17 


0.14 


Rm 


0.194 


0.197 


0.199 


0.201 



Table 1: The PQCD estimations about the hadronic matrix element ( K \ s'~^^h \ B ) and 
the ratio of ( r;' | S75S | 0) ( K~ \ sb \ B^) to ( 77' | 57^755 | 0) ( K~ \ s'y^b \ B^). 
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BTZ(B~ n'K~\ 


Merh Cw) 
contribution 


IVTppIi riv^-l-rroTivpTitioTiall 
contribution 


Mech. (iv) + (conventional) +Mech. (iii) 


f^, = lOMeV 


r, = -15Mey 


= -30MeV 


€B = 0.05 


2.65 X 10"^ 


4.67 X 10-5 


4.31 X 10-5 


5.45 X 10-5 


6.51 X 10-5 


(11 = {).{)() 


2.0 X K)-"' 


.S.<sr) X K)-"' 


.S.49 X K)-"' 


i.ir.i X 


.1.70 X K)-" 


CB = 0.07 


1.57 X 10-5 


3.32 X 10-5 


2.96 X 10-5 


4.10 X 10-5 


5.17 X 10-5 


CB = 0.08 


1.30 X 10-5 


2.96 X 10-5 


2.60 X 10-5 


3.74 X 10-5 


4.80 X 10-5 



Table 2: Numerical results of BTZ{B rj'K ) in different cases. For conventional mecha- 
nism(tree+penguin), Bn{B- rj'R-) = 0.9 x 10"^ 



Branching 

riatio 


Mech.(iv) 

(B = ■()(')({). i)7) 


€3 = 0.06(0.07) 

f;,. = -i'>M<v 


Ref[3] 


Ref[10] 


Ref[3]+ 

Mocli.(iv) 


Exp.(10-5) [1] 


B- rj'K- 


2.0(1.57) X 10-5 


4.63(4.10) X 10-5 


5.69 X 10-5 


6.1 X 10-5 


8.31 X 10-5 


6.5tij±0.9 


B^ ^ r]'K° 


2.0(1.57) X 10-5 


4.03(3.66) X 10-5 


5.19 X 10-5 


5.69 X 10-5 


6.24 X 10-5 


AJtl'l ± 0-9 


B- rj'K*- 


1.52(1.12) X 10-5 


3.18(2.61) X 10-5 


3.24 X 10-^ 


1.15 X 10-5 


7.41 X 10-*^ 


< 13 


B° T]'K*" 


1.52(1.12) X 10-5 


3.33(2.74) X 10-5 


2.96 X 10-^ 




7.29 X 10-6 


< 3.9 


B- r]K- 


3.46(2.77) X 10-*^ 


5.31(4.52) X 10-6 


8.68 X 10-^ 


5.13 X 10-6 


5.29 X 10-6 


< 1.4 




3.46(2.77) X 10-6 


5.02(3.90) X 10-6 


3.28 X 10-^ 


4.72 X 10-6 


4.80 X 10-6 


< 3.3 


B- ^ rjK*- 


2.76(2.05) X 10-6 


3.44(4.20) X 10-6 


3.70 X 10-6 


1.31 X 10-5 


9.73 X 10-^ 


< 3.0 


B° r)K*° 


2.76(2.05) X 10-6 


2.40(3.09) X 10-6 


2.44 X 10-6 




4.21 X 10-"^ 


< 3.0 



Table 3: Numerical results for B decay channels of interest and compared with predictions 
of other references and experimental results. 
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Figure 2: Leading twist diagrams for {K \ 57^6 \ B ) in PQCD. 
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Figure 3: Leading twist diagrams for {K \ sb \ B ) in PQCD. 
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Figure 4: Diagram for the two gluons fusion mechanism. 
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